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Abstract 

A search for a Higgs boson decaying into invisible particles is performed using the 
data collected at LEP by the L3 experiment at centre-of-mass energies of 183 GeV 
and 189 GeV. The integrated luminosities are respectively 55.3 pb~^ and 176.4 pb~^. 
The observed candidates are consistent with the expectations from Standard Model 
processes. In the hypothesis that the production cross section of this Higgs boson 
equals the Standard Model one and the branching ratio into invisible particles is 
100%, a lower mass limit of 89.2 GeV is set at 95% confidence level. 
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1 Introduction 



In some extensions of the Standard Model the Higgs boson can decay into stable weakly inter- 
acting particles, thus yielding invisible final states |1J. For example the minimal supersymmetric 
extension of the Standard Model predicts that the Higgs boson can decay into a pair of invisible 
neutralinos. 

A search is performed for a Higgs boson produced through the Higgs-strahlung process, 
e"'"e~ — > Z* — > hZ. The Z boson decays into fermion pairs yielding two different investigated 
topologies: two acoplanar jets plus missing energy, corresponding to the Z boson hadronic 
decays, and two acoplanar charged leptons plus missing energy, corresponding to decays of 
the Z boson into electrons or muons. Data collected by the L3 experiment |2j at LEP centre- 
of-mass energies of ^/s = 183 GeV and 189 GeV are analysed. The corresponding integrated 
luminosities are respectively 55.3 pb~^ and 176.4 pb~^. Results at lower centre-of-mass energies 
have been reported by L3 3] and by the other LEP experiments j3]. 

2 Event simulation 

To determine the signal efficiency, samples of Higgs boson events are generated using the 
PYTHIA Monte Carlo program [S] for masses between 55 GeV and 100 GeV. 

For the background studies the following Monte Carlo programs are used: PYTHIA (e"'"e~ 
qq(7), e+e- ^ Z/7*Z/7* and e+e" ^ Ze+e"), KORALW (e+e" W+W), KORALZ [Zj 
(e+e- ^ yU+/i-(7), e+e- ^ t+t-{j)), PHOJET jBl (e+e" ^ e+e-qq), DIAG36 ^ (e+e" ^ 
e+e" £+£-), BHWIDE ^ (e+e" ^ e+e" (7)), and EXCALIBUR ^ for the other four-fermion 
final states. For each centre-of-mass energy, the number of simulated background events corre- 
sponds to at least 50 times the number of expected events except for the two-photon interactions 
(e+e" — >■ e+e~ff) and Bhabha scattering (e+e~ — > e+e~) for which twice and seven times the 
collected luminosity are simulated, respectively. 

The L3 detector response is simulated using the GEANT 3.15 program jT2], which takes 
into account the effects of energy loss, multiple scattering and showering in the detector. The 
GHEISHA program is used to simulate hadronic interactions in the detector. Small time- 
dependent inefficiencies of the different subdetectors are taken into account in the simulation 
procedure. 

3 Search in the hadronic channel 

A cut-based analysis is used to select events in the hadronic channel. After a common prese- 
lection two sequential selections are separately optimised for light (below 80 GeV) and heavy 
(above 80 GeV) Higgs boson masses. Unless otherwise stated, the events are constrained to 
have two jets using the DURHAM algorithm p3|- 

3.1 Preselection 

High-multiplicity hadronic events are selected at ^/s = 189 GeV. Events coming from QCD 
processes and hadronic decays of W and Z boson pairs are rejected by requiring a missing 
momentum larger than 10 GeV. The absolute values of the cosine of the polar angle of the 
jets and of the missing momentum vector have to be less than 0.9, to reject events with a 
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high-energy initial-state radiation photon emitted close to the beam axis. In addition, events 
with large energy depositions in the forward calorimeters are vetoed to reduce the background 
contribution from the Ze"'"e~ and qq(7) processes and from residual two-photon interactions. 
Coplanar events are rejected to further suppress these last processes. 

Events with energetic and isolated charged leptons are removed to decrease the contami- 
nation from semi-leptonic decays of W boson pairs. This cut is designed to keep signal events 
with semi-leptonic decays of b or c hadrons produced in Z decays. 

The larger of the jet masses is required to be in the range from 6 GeV up to 50 GeV and 
the lower one greater than 4 GeV. The upper mass limit further removes some semi-leptonic 
W boson pair decays and the lower bounds reject two-photon interactions with tau leptons in 
the final state. 

Figure n shows the comparison between data and Monte Carlo expectations for the distri- 
bution of the visible mass, M^is of the preselected events. 

3.2 Heavy Higgs boson selection 

In addition to the preselection described above, the correlation between the visible and the 
missing mass, Mmis, is used to select heavy Higgs boson candidates in the y/s = 189 GeV data 
sample. We define the variable R = (Mvis + Minis)/(Mvis — Mmis) and we require R < —3.5 or 
R > 7, since the signal has a broad R distribution while for the background R is close to zero. 

A heavy Higgs boson is characterised by relatively low momentum, hence the missing mo- 
mentum of the event should not exceed 40 GeV. The background due to the qq(7) and two- 
photon interaction processes is suppressed by rejecting coUinear events and by requiring a large 
value of the event thrust, together with a moderate value of the sum of the inter-jet angles, 
Oi23, when the events are constrained to have three jets. In order to reject the residual con- 
tributions from the W pair and the single W processes, an isolation criterion on the missing 
momentum vector is applied. In addition, the upper cut on the maximum value of the jet 
masses is tightened to 30 GeV while the minimum has to be less than 20 GeV. 

The recoiling mass, M^'^'^, is calculated by constraining the visible mass to the Z boson mass 
and imposing energy-momentum conservation its distribution is plotted in Figure Efa) for 
the data and the background. With this kinematical constraint, the recoil mass resolution is 
3.5 GeV in the hypothesis of a Higgs boson mass of 90 GeV. 

After applying the selection described above the dominant process in the remaining back- 
ground is the Z boson pair production. 

3.3 Light Higgs boson selection 

The production of a light invisible Higgs boson at ^/s = 189 GeV is characterised by three 
main features, exploited by the following selection criteria: mass of the hadronic system close 
to the Z mass |Mvis — Mz| < 20 GeV, at least 40% but not more than 60% of the centre-of- 
mass energy visible in the detector and missing momentum in the window from 30 GeV up 
to 55 GeV. This last requirement reduces part of the background arising from Z boson pair 
production and two-photon interactions, the latter being further suppressed by an upper cut 
on 0123. Events from qq(7) are rejected by requiring a large value of the event thrust and 
the longitudinal momentum imbalance to be less than 40% of the visible energy. The residual 
contribution from W pair production is reduced by a cut on the threshold j/23 at which the 
DURHAM algorithm resolves the event into three jets from a two-jet topology. 
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The distribution of M^'^'^ for events selected in the data and the Monte Carlo samples is 
displayed in Figure El^b). 

The selection of the hadronic channel at ^/s = 183 GeV is similar to this light Higgs boson 
selection. A cut on the transverse momentum imbalance replaces the cut on the 1/23 parameter 
and the values of the selection requirements reflect the different centre-of-mass energies. The 
final M^^^^ spectrum for data and Monte Carlo is shown in Figure |2l^c). Tabled] summarises the 
yields of all the selections described above. 

After applying the two selections described above the dominant process in the remaining 
background is the W"'"W~ production. 
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Table 1: Number of events expected from Standard Model processes compared to 
the number of data events selected by the hadronic selections. The signal efficiencies 
(e) for several Higgs boson masses Mh are also shown. 



4 Search in the leptonic channels 

The search for an invisibly decaying Higgs boson produced in association with a Z boson decay- 
ing into leptons is designed to be almost independent of the Higgs boson mass in the investigated 
range. Low multiplicity events with a pair of high energy muons or electrons are selected. These 
are separated from fermion pair production events by requiring large acoplanarity and visible 
energy between 5% and 70% of the centre-of-mass energy. The lepton energy has to be less than 
90% of the beam energy to further reject Bhabha scattering events. Two-photon interactions 
are suppressed by requiring the lepton pair invariant mass to be larger than 30 GeV and low 
energy depositions in the forward calorimeters. Events with muons should have at least one 
scintillator in time with the beam crossing in order to remove cosmic-ray background. The 
yield of this preselection is presented in Table El while Figure El displays the spectra of the 
lepton pair invariant mass, M^, for data and Standard Model Monte Carlo events. 

Residual events due to the radiative return to the Z resonance where the photon remains 
undetected in the beam pipe are rejected by requiring the missing momentum to point away 
from the beam axis. Tau pair production can yield acoplanar lepton pairs that satisfy the 
selection criteria described above. In the hypothesis that the lepton pair originates from a 
single particle, we require the cosine of the most energetic lepton emission angle 6* in the Z 
boson rest frame not to exceed 0.95. 

The contribution from two-photon interactions is eliminated by tightening the cut on the 
lepton invariant masses, 70 GeV < M« < 110 GeV; this is also effective against a significant 
portion of the fully leptonic decays of W bosons. Final states with an electron or muon pair 
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Table 2: Number of events observed and expected from Standard Model processes 
at -^i = 189 GeV after the preselections and the final selections. Signal efficiencies 
(e) for different Higgs mass hypotheses are also shown. The background to the final 
selected sample is composed of one third Z boson-pair events and two thirds W 
boson-pair events. 



and two neutrinos, produced by Z boson pairs, constitute an irreducible background but their 
cross section is relatively low. 

The visible energy E'vis, cos M^^ and the velocity /5 of the dilepton system are combined 
into a single likelihood variable G, defined as: 

G = Y.^\0^{Pl{x))-\0^[P],{x)). 

The index % runs over the four variables and Ps{x) and Pb{x) are the probability densities 
for the i-th variable to have a value x in the signal or background hypotheses, respectively. 
These densities are calculated for each event by interpolating between the two signal Monte 
Carlo samples whose generated Higgs masses are closer to the event missing mass which is 
taken as the Higgs boson mass hypothesis. The Z boson pair background is not included in 
this calculation. Figure ID shows distributions of G for the data and the expected Monte Carlo 
background and signal for a Higgs boson mass of 95 GeV. 

The number of selected events and the signal efficiency after the optimization ^3] of a cut 
on G are reported in Table |21 The observed resolution on the missing mass is 1.1 GeV in the 
electron channel, and 5.1 GeV in the muon channel for a Higgs boson mass of 90 GeV. 

A cut-based analysis is developed for the ^/s = 183 GeV data sample making use of the 
following selection criteria: 30 GeV < E^-,, < 120 GeV, cos 9* < 0.95, 80 GeV < M« < 
100 GeV and 0.05 < (3 < 0.55. In the electron channel, the signal efficiency is 45% and 4 events 
are observed for 1.4 expected background events. In the muon channel, the signal efficiency is 
28% and no events are observed while 1.7 background events are expected. Figure El displays the 
missing mass distributions for the search in the leptonic channel for the combined ^/s = 183 GeV 
and y/s = 189 GeV samples. 

5 Systematic uncertainties 

Two sources of systematic uncertainties, summarised in Table El can affect the results. The 
first is the limited amount of Monte Carlo statistics, which gives the systematic errors on 
the signal and background efficiencies listed as "MC Stat." in Table El The second is the 
quality of the Monte Carlo description of the background processes. This is studied using data 
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and Monte Carlo samples containing essentially W^W~ and ZZ background events. These 
samples contain about 1100 events for the hadronic channel and 500 for the leptonic ones. The 
data distributions in these new samples of each selection variable i, except the likelihood G, 
are compared with those of the Monte Carlo, determining their systematic shifts Si and the 
corresponding statistical errors Uj. All the selection cuts are then shifted by Sj ± cxj, where the 
sign of (jj is chosen so as to obtain the lowest efficiency for the single cut on the variable i. The 
difference between the efficiency of the selection using the shifted cuts and that of the nominal 
one is taken as the systematic uncertainty. These errors are summarised as "Syst." in Table IHl 
and are summed in quadrature with the Monte Carlo statistical uncertainties to obtain the 
total systematic uncertainty, listed as "Total" in Table El 
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Table 3: Relative systematic uncertainties in percent on the signal and background 
efficiencies for each analysis. 



6 Results 

No indication of the production of a Higgs boson with invisible decays is found. As both the 
production cross section and the branching ratios are model dependent, it is useful to introduce 
the ratio Ri^v = BR(h — > invisible particles) x (T(e+e~ hZ)/a{e'^e~ — >• HsmZ), where Hsm 
is the Standard Model Higgs boson. A limit on i?inv is calculated (TH] as a function of the 
Higgs boson mass making use of the mass distributions presented in Figures |21 and In 
the determination of the limit the Standard Model Higgs boson cross section as given by the 
HZHA generator ^7] is used and the signal and background efficiencies are lowered by their 
systematic uncertainties. Results obtained at lower energies |S] are included. Figure IHl shows 
the 95% confidence level (CL) upper limit on Rmv as a function of the Higgs mass Mh. For the 
value of i?inv = 1 the 95% CL lower limit on the Higgs boson mass is: 

Mh > 89.2 GeV. 

The expected lower limit is 92.6 GeV. 
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Figure 1: The distribution of the visible mass after the hadronic preselection at 
= 189 GeV. 
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Figure 2: The distribution of the recoil mass (a) after the heavy Higgs boson 
selection at y/s = 189 GeV, (b) after the light Higgs boson selection at ^/s = 
189 GeV and (c) after the final selection at i/i = 183 GeV. 



11 




M,, (GeV) 



8' 



4 DATA 

D Two-photon process 
D e"^e" \x^\L (y) 
a WW, ZZ, Zee 



(b) 




M,, (GeV) 



120 140 



Figure 3: Distribution of (a) the dielectron and (b) the dimuon invariant mass at 
y/s = 189 GeV after the preselection is apphed. 
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Figure 4: Distributions of the final likelihood variable, G, (a) for the electron and 
(b) for the muon selections at ^/s = 189 GeV, for data and the expected background. 
A possible Higgs signal (Mh = 95 GeV) with an arbitrary cross section is also shown. 
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Figure 5: The missing mass distributions (a) in the electron channel and (b) in the 
muon channel for the combined ^/s = 183 GeV and 189 GeV data samples. 
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Figure 6: Observed and expected upper limits on the ratio of the invisibly-decaying 
Higgs boson cross section to that of the Standard Model Higgs boson, as function 
of the Higgs boson mass. The shaded area is excluded at least at 95% confidence 
level. 



15 



